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ABSTRACT 

We present maps of ^^CO J = 2 — 1 emission covering the entire star- forming disks of 16 nearby 
dwarf galaxies observed by the IRAM HERACLES survey. The data have 13'' angular resolution, 
~ 250 pc at our average distance of D = 4 Mpc, and sample the galaxies by 10 — 1000 resolution 
elements. We apply stacking techniques to perform the first sensitive search for CO emission in 
dwarf galaxies outside the Local Group ranging from individual lines-of-sight, stacking over IR-bright 
regions of embedded star formation, and stacking over the entire galaxy. We detect 5 galaxies in CO 
with total CO luminosities of Leo 2-1 = 3 — 28 x 10^ K km s~^ pc^. The other 11 galaxies remain 
undetected in CO even in the stacked images and have Leo 2-1 ^ 0.4 — 8 x 10^ K km s~^ pc^. We 
combine our sample of dwarf galaxies with a large sample of spiral galaxies from the literature to 
study scaling relations of Leo with Mb and metallicity. We find that dwarf galaxies with metallicities 
of Z ?^ 1/2 — 1/10 Zq have Leo of 2 — 4 orders of magnitude smaller than massive spiral galaxies 
and that their Leo per unit Lb is 1 — 2 orders of magnitude smaller. A comparison with tracers of 
star formation (FUV and 24/im) shows that Leo per unit SFR is 1 — 2 orders of magnitude smaller 
in dwarf galaxies. One possible interpretation is that dwarf galaxies form stars much more efficiently, 
we argue that the low Leo/SFR ratio is due to the fact that the CO-to-H2 conversion factor, aco, 
changes significantly in low metallicity environments. Assuming that a constant H2 depletion time 
of Tdep = 1.8 Gyr holds in dwarf galaxies (as found for a large sample of nearby spirals) implies 
aeo values for dwarf galaxies with Z ^ 1/2 — 1/10 Zq that are more than one order of magnitude 
higher than those found in solar metallicity spiral galaxies. Such a significant increase of aeo at 
low metallicity is consistent with previous studies, in particular those of Local Group dwarf galaxies 
which model dust emission to constrain H2 masses. Even though it is difficult to parameterize the 
dependence of aeo on metallicity given the currently available data the results suggest that CO is 
increasingly difficult to detect at lower metallicities. This has direct consequences for the detect ability 
of star-forming galaxies at high redshift which presumably have on average sub-solar metallicity. 

Subject headings: galaxies: ISM — ISM: molecules — radio lines: galaxies 



1. INTRODUCTION 

Robust knowledge of the molecular (H2) gas dis- 
tribution is indispensable to understand star forma- 
tion in galaxies. Observations in the Milky Way 
and nearby galaxies suggest that stars form in clouds 
consisting predomin antly of H2 (|Lada fc Ladal 120031 : 
iFukui fc Kawamiiral[2010.) . Because H2 is almost impos- 
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sible to observe directly under typical conditions of the 
cold interstellar medium (ISM), its abundance and distri- 
bution has to be inferred using indirect methods. Obser- 
vations of low rotational lines of carbon monoxide (CO) 
have been the standard method to do so as CO is the 
second most abundant molecule and easily excited in the 
cold ISM. Over the last decades, of the order of a hundred 
galaxies in the local Universe have been successfully de- 
tected in CO. Over the last years CO has been detected 
throug:hout the Universe out t o cosmological distances 
(jSolomon fc Vanden Boutll2005f ). These CO observations 
have greatly enhanced our knowledge of H2 in galaxies, 
the phase balance of the ISM, and its interplay with star 
formation. 

Despite great advances in studying H2 in massive star- 
forming galaxies, our knowledge of H2 in star- forming 
dwarf galaxies remains poor. The CO emission in these 
systems has proven to be extremely faint and most stud- 
ies targeting metal-poor dwarf galaxies have resulted 
in non-detections. For sensitivity reasons, surveys of 
dwarf galaxies have tended to target only a few sys- 
tems and used mostly sing:le point ing:s (llsrael et al.[[T995 



igs^ 

'Young et aLi'iggS"; 'Tavlor et aLi'iggS"; 'Barone et al."2^; 
Boker et al. 2003; Sauty et al. 2003; Albrecht et al. 2004; 
ILerov et al.ll2005l) . These data are very heterogeneous as 
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they target different CO transitions, cover different re- 
gions, and have different beam sizes, sensitivities, and 
beam fihing factors. Thus, conclusive results for ba- 
sic quantities such as the total CO luminosity of dwarf 
galaxies have not been reached and comparison to other 
observables have been complicated by these systematic 
effects. 

CO observations are currently — and will remain — 
our most accessible tracer of cold H2 in the local and 
distant Universe. It is thus important to obtain pro- 
found understanding of the connection between CO and 
H2 in different environments. Inside individual molecular 
clouds, this dependence has proven to be h ighly compli- 
cated and influenced by many factors fe.g. JShettv et al.l 
l2Qllallbl ). Many of these dependencies average out on 
scales larger than individual clouds, however metallic- 
ity will not. Metallicity may thus be the single most 
important factor determining the CO/H2 ratio on large 
scales. This makes a robust calibration of the CO/H2 
ratio as function of metallicity a viable proposition for 
those observational studies that use CO as a tracer of 
H2. The need becomes more pressing as observations 
start probing the CO content of galaxies in the distant 
Universe where most stars presumably formed in envi- 
ronments with sub-solar metallicity. 

To understand the environmental dependencies of the 
CO/H2 ratio requires good knowledge of the CO content 
of all types of galaxies, even in those where we worry 
that CO may not trace H2 in the same way as it does 
in massive spiral galaxies. This makes sensitive, wide- 
field CO maps of dwarf ga laxies an important under- 
takin^. The HER ACLE^ survey (partly pubhshed in 
iLerov et"aI]|2QQ9l ) has obtained such CO observations of 
a large set of nearby star-forming galaxies ranging from 
massive spirals down to low-mass, low- metallicity dwarfs. 
In conjunction with an extensive set of multi- wavelength 
data, this survey has already led to a vast improvement 
of our knowledge of the relation between Hi, CO, H2, 
and star formation. 

In this paper, we present sensitive measurements of 
CO emission of 16 nearby low-mass, low- metallicity star- 
forming dwarf galaxies from the HERACLES survey us- 
ing stacking techniques. We use these data to study the 
relation between CO emission and other galaxy parame- 
ters, especially star formation rate (SFR) and H2 mass. 
Then we analyze the CO/H2 ratio as function of metal- 
licity. In Section [2] we introduce our multi- wavelength 
data and summarize their basic properties. In Section [3] 
we conduct a sensitive search for CO emission for indi- 
vidual lines-of-sight, IR-bright regions, and entire galax- 
ies. In Section |4] we compare these CO measurements to 
other galaxy parameters and compare the relationships 
found for dwarf galaxies to those of massive spiral galax- 
ies. In Section [5] we study the metallicity dependence 
of the CO/H2 ratio. We use observed SFRs to infer H2 
masses and thus constrain CO/H2, then we compare our 
results to results derived from other methods. In Sec- 
tion [6] we summarize our findings. 

2. DATA 

We study 16 nearby low-mass star- forming galaxies (for 
simplicity just called "dwarfs" throughout the paper) 

^^ |http : //www . cv . nrao . edu/~ aleroy/HERACLES/Overview . html | 



from the HERACLE^ survey (|Lerov et al.l [20091 ): see 
Figure [1] for an outline of the area surveyed for each tar- 
get. These data have been largely neglected in previous 
work on the HERACLES sample as their CO emission 
has rarely been rob ustly detected in t he pixel-based stud- 
ies of Bigie l et al.l (2008, 2011) a nd ILerov et all (I2008D 
or the radial stacking analysis of iSchruba et al.l (|2011[ ). 
The only galaxies for which data were provided are Ho I, 
Ho II, DDO 154 IC 2574, NGC 2976, and NGC 4214 
(Leroy et ajj l2009f ) . Table [1] lists our sample of dwarf 
galaxies along with adopted distances, inclination, po- 
sition angle, optical radius, metallicity, 5-band opti- 
cal magnitude. Hi mass, and total sta r formation rate 
(SFR). These values are taken from W alter et al.l (I2008D 
where possible and from LEDA (;Prugnlei fc Heraudeaul 
Il998f ) and NED elsewhere. 

2.1. CO Data 

We take CO data from the HERACLE^ survey which 
mapped the ^^CO J = 2 ^ 1 emission line i n 48 nearby 
galax ies using the IRAM 30m telescope ( Lerov et al.l 
l2009f ) . The observations are designed to cover large parts 
of the galaxies and extend to 1 — 1.5 times the optical 
radius, i?25, for large spirals and up to 2 — 3 i?25 for small 
galaxies. The final data cubes have an angular resolution 
(FWHM) of 13^' and a spectral resolution (channel sepa- 
ration) of 2.6 km s~^. The noise level is 20 — 30 mK per 
resolution element and per channel. 

Whenever possible we compare our observed CO (2-1) 
data to literature measurements. To do that, we convert 
them to CO (1-0) intensities assuming a constant line ra- 
tio, R21 = Ico 2-1 / ^co 1-0 = 0.7. We choose this con- 
stant line ratio to achieve consi stency with Bigiel et al.l 
(|2011f ) and ISchruba et al.l (j2011f ). This is the average ra- 
tio found for all HERACLES galaxies (E. Rosolowsky et 
al., in preparation). From this data no significant varia- 
tions of R21 with metallicity are evident over the range 
12 + log 0/H ^ 8.6 — 8.9 on scales of 1 kpc. In low metal- 
licity environments such as the Magellanic clouds values 
of i^2i ^ 1 -0 — 1.5 are frequently found (jBolatto et al.l 
I2000ri2003l : llsrael et al.ll2003Ll2005f ). This can lead to an 
over-prediction of the true CO (1-0) values by a factor of 
< 2. However, later we will see that this potential bias 
is too small to change our conclusions. 

We will discuss the CO-to-H2 conversion factor 
extensively in Section [5] but note that a typical 
Galactic CO(l-0)-to-H2 conversion factor is Xqo = 
2.0 X lO^Q cm-^ f K km s-^)-^ (ISta^g fc Mattoxl [19961 : 
iDame et al.l 120011 : lAbdo et al.ll2010f ) which translates to 
<^co = 4.38 M0 pc~^ (K km s~-^)~^ when including a 
factor of 1.36 to account for heavy elements. 

2.2. HI Data 

W e draw H I data mos tly from the VLA THINGS sur- 
vey (jWalter et al.l[2QQ8h . The Hi data for NGC 4236, 
NGC 4626, and NGC 5474 are from the VLA programs 
AL 731, AL 735, both led by one of us (P.I. Leroy), and 
from the archive. The Hi data for DDO 165 is from the 
LITTLE THINGS survey (D. Hunter et al., in prepa- 
ration). The angular resolution of the data referred to 
above is ~ 10 — 20'^ the velocity resolution is 2.6 — 5.2 
km s-i for the THINGS and LITTLE THINGS data, 
and 5.2 — 10 km s~^ for the other data. The sensitiv- 
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Figure 1. Coverage of CO Observations. For each HERACLES dwarf galaxy we show the coverage of our CO data (gray dashed hne), the 
Hi surface density (grayscale) at linear scale between — 40 Mq pc~^, the 24/im intensity (red contour) at 0.2 MJy sr"-*^, the FUV intensity 
(blue contour) at 0.01 MJy sr"-*^, and a galactocentric radius R = R25 (white contour). We determine the CO intensity at each sampling 
point (dot), and stack the data for the "entire" galaxy (black contour) and for IR-bright regions (red contour), see text for definition of 
stacking regions. 
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Figure 1. (Continued) 
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Table 1 

Properties of Galaxy Sample 



Name 


Alt. Name 


D Incl. 


P.A. 


R25 


Metal. ^ 


Mb^ 


Mhi^ 


SFRd 






(Mpc) n 


n 


n 


12+logO/H 


(mag) 


(108 Mo) 


(Mo yr-i) 


M81 DwA 




3.6 23 


49 


0.64 


7.50 


-11.4 


0.12 


0.0005 


M81 DwB 


UGC 5423 


5.3 44 


321 


0.56 


8.02 


-13.8 


0.25 


0.0023 


DDO 053 


UGC 4459 


3.6 31 


132 


0.39 


7.80 


-13.9 


0.60 


0.0035 


DDO 154 


UGC 8024 


4.3 66 


230 


0.98 


7.78 


-15.4 


3.58 


0.0056 


DDO 165 


UGC 8201 


4.6 51 


90 


1.66 


7.84 


-14.1 


6.33 


0.0100 


HO I 


UGC 5139 


3.8 12 


50 


1.65 


7.83 


-16.8 


1.39 


0.0100 


HO H 


UGC 4305 


3.4 41 


177 


3.76 


7.93 


-12.5 


5.95 


0.0455 


IC 2574 


UGC 5666 


4.0 53 


56 


6.41 


8.05 


-17.2 


14.80 


0.0718 


NGC 2366 


UGC 3851 


3.4 64 


40 


2.20 


7.96 


-16.2 


6.49 


0.0605 


NGC 2403 


UGC 3918 


3.2 63 


124 


7.87 


8.57 


-18.6 


25.80 


0.4140 


NGC 2976 


UGC 5221 


3.6 65 


335 


3.60 


8.67 


-16.5 


1.36 


0.0895 


NGC 3077 


UGC 5398 


3.8 46 


45 


2.70 


8.64 


-17.3 


8.81 


0.0838 


NGC 4214 


UGC 7278 


2.9 44 


65 


3.40 


8.25 


-17.1 


4.08 


0.1208 


NGC 4236 


UGC 7306 


4.4 75 


162 


11.99 


8.46 


-18.1 


34.60 


0.1409 


NGC 4625 


UGC 7861 


9.5 47 


330 


0.69 


8.70 


-17.0 


11.80 


0.0716 


NGC 5474 


UGC 9013 


6.8 50 


85 


1.20 


8.57 


-17.3 


15.50 


0.1069 


^ OxvEcen abundance from IMoustakas et al 


.1 (|2010l) 












^ B-band maj 


mitude from HERACLES. 














^ Mmt fromlWalter et al.l (| 


20081). 















d SFR(FUV+24) from this work. 

ity of these data is sufficiently high to never hmit our 
analysis. 

2.3. Star Formation Tracers 

We estimate the star formation rate (SFR) using a 
combination of FUV and 24/im emission follo wing the ap- 
proach introduced in Bigiel et al. (2008) and lLerov et al.l 
20081 ). The SFR surface density is given by User 
[Mo yr-i kpc-2] = 0.081 (/fuv + 0.04 /24;.m) [MJy sr-^] 
X cos i. The FUV dat a are taken froni the GALEX 
Nearby Galaxy Survey (|Gil de Paz et al.1 120071 ) or al- 
ternatively from the NASA Multimission Archive at 
STScI. They cover a wavelength range of 1350 — 1750 A, 
have angular resolution ~ 4.5'', and sufficient sensitiv- 
ity to determine FUV intensities with high signal-to- 
noise throughout the star-forming disk. The IR data are 
taken from the Spitzer SINGS (Ken nicutt et al.' "20031) 
and Local Volume Legacy (LVL) surveys (jPale et al. 
I2OO9). These data have ~ 6'' resolution; their sensitiv- 
ity is sufficient to detect 24/im emission in most of our 
galaxies except the lowest mass and lowest metallicity 
dwarf galaxies. We apply some processing to the FUV 
and 24/im maps (i.e., mask foreground stars and flatten 
background) as described in ILerov et al.l (|2012f ). 

2.4. Metallicities 

Ga s phase oxygen ab u ndanc es (metallicities) are taken 
from IMoustakas et al.l ()2010l ). For galaxy-integrated 
data we use the average of their "characteristic" metal- 
licities (their Table 9) derived from a theoretical calibra- 
tion (KK04 values) and an empirical calibration (PT05 
values). In s ome plots we also sh ow the radial stacking 
results from Schruba et al.l (120111) for which we us e the 
metallicity gradients from IMoustakas et al.l ()2010L Ta- 
ble 8) again averaging the two calibrations. The "charac- 
teristic" metallicity of a galaxy corresponds to the value 
of the metalHcity gradient at radius R = 0.4 R25. 

2.5. Sampling 



We convolve all our data to a common resolution of 
13'' (limited by the CO data; some Hi data have coarser 
beam sizes and we include them on their native resolution 
assuming to ffist order homogeneous Hi distribution) 
and sample them on a hexagonally packed grid spaced 
by half a beam size (6.5"). For each line-of-sight we col- 
lect observed intensities of CO, H I, FUV, and 24/im, and 
determined local gas masses and SFRs. We also store 
the H I mean velocity, the original CO spectrum, and the 
galactocentric radius. Figure [T] shows for each galaxy in 
our sample the Hi distribution as grayscale, the extent 
of FUV and 24/im emission indicated by a single con- 
tour, the CO map coverage, and our sampling grid as 
dots. See Figure |2] for integrated CO intensity maps for 
a subset of our galaxies. 

2.6. Literature Sample 

Throughout the paper we will compare our measure- 
ments for dwarf galaxies to a larger sample of nearby 
galaxies. T his sample is tak e n from the literature com- 
pilation of iKrumholz et al.l ()2011l ) which includes the 
more massive HERACLES galaxies and some additional 
Local Group and nearby galaxies. Table [2] lists their 
names together with adopted distances, metallicities, B- 
band magnitudes, total CO (1-0) luminosity, total SFR, 
and references to the original literature. The compila- 
tion aims at maximizing homoge neity of used dat a and 
methodology. We supplement the IKrumholz et aTl com- 
pilation by adding absolute 5-band magnitudes adjusted 
to our adopted distances. We also update the total 
CO luminosities using the most recent HERACLES data 
(converted to CO (1-0) luminosities) and SFRs derived 
from combinin g FUV and 24/im maps. We use the metal- 
licities listed in IKrumholz et ahl . these have been derived 
following the above described methodology. 

3. CO EMISSION IN HERACLES DWARF GALAXIES 

To derive meaningful constraints on CO content, we 
search for CO emission on three different spatial scales: 
individual lines-of-sight, stacked over the entire galaxy 
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Figure 2. CO Intensity Maps. Integrated CO(2-l) intensity maps (including all channels with velocities in the range of ±50 km s ^ 
of the local Hi mean velocity) for a subsample of the HERACLES dwarf galaxies at a linear grayscale from —0.3 to 0.7 K km s"-*^. The 
galactocentric radius R = R25 is shown as white contour, the region selected to determine the integrated CO intensity for the entire galaxy 
as black contour. The remaining HERACLES dwarf galaxies not shown here are non-detections at full resolution. 
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Table 2 

Properties of Literature Galaxy Sample 



Name 



D 

(Mpc) 



Ref.^ 



Metal. 
12+logO/H 



Ref/ 



Mb 
(mag) 



Ref.^ 



log Leo 1-0 
(K km s-i pc2) 



Ref.^ 



log SFR 

(M0 yr-i) 



Ref.^ 



SMC 


0.06 


Lll 


8.00 


D84; MAIO 


-16.2 


Lll 


5.20 


M06 


-1.30 


W04 


LMC 


0.05 


Lll 


8.30 


D84; MAIO 


-17.6 


Lll 


6.50 


F08 


-0.70 


H09 


IC 10 


0.95 


HOI 


8.20 


L79; L03 


-16.5 


HOI 


6.30 


L06 


-1.03 


L06 


M 33 


0.84 


G04 


8.30 


R08 


-18.9 


NED 


7.60 


H04 


0.00 


H04 


I ZW 18 


14.00 


104 


7.22 


T05 


-14.7 


G03 


< 0.10 


L07 


-1.00 


L07 


II ZW 40 


9.20 


CIO 


8.10 


E08; CIO 


-17.9 


NED 


6.20 


T98 


-0.19 


CIO 


NGC 0628 


7.30 


W08 


8.69 


MOlO 


-20.0 


W08 


8.45 


HERA 


-0.08 


HERA 


NGC 0925 


9.20 


W08 


8.52 


MOlO 


-20.0 


W08 


7.47 


HERA 


-0.24 


HERA 


NGC 1482 


22.00 


CIO 


8.53 


MOlO 


-18.8 


NED 


8.80 


Y95 


0.53 


CIO 


NGC 1569 


3.36 


G08 


8.10 


M97 


-18.1 


NED 


5.55 


T98 


-0.40 


Pll 


NGC 2146 


12.80 


W08 


8.70 


E08; CIO 


-20.6 


W08 


9.06 


HERA 


0.93 


CIO 


NGC 2537 


6.90 


Lll 


8.40 


MAIO 


-16.4 


Lll 


5.50 


T98 


-1.05 


CIO 


NGC 2782 


40.00 


CIO 


8.60 


E08; CIO 


-20.9 


NED 


9.00 


Y95 


0.72 


CIO 


NGC 2798 


24.70 


W08 


8.69 


MOlO 


-19.4 


W08 


8.73 


HERA 


0.49 


CIO 


NGC 2841 


14.10 


W08 


8.88 


MOlO 


-21.2 


W08 


8.34 


HERA 


-0.10 


HERA 


NGC 2903 


8.90 


W08 


8.90 


MAIO 


-20.1 


Lll 


8.82 


HERA 


0.32 


HERA 


NGC 3034 


3.60 


W08 


8.82 


MOlO 


-18.5 


Lll 


8.94 


HERA 


0.90 


CIO 


NGC 3079 


21.80 


CIO 


8.60 


E08; CIO 


-21.7 


NED 


9.40 


Y95 


0.50 


CIO 


NGC 3184 


11.10 


W08 


8.83 


MOlO 


-19.9 


W08 


8.56 


HERA 


-0.01 


HERA 


NGC 3198 


13.80 


W08 


8.62 


MOlO 


-20.7 


W08 


8.15 


HERA 


-0.01 


HERA 


NGC 3310 


21.30 


CIO 


8.20 


E08; CIO 


-20.5 


NED 


8.20 


Y95 


0.92 


CIO 


NGC 3351 


10.10 


W08 


8.90 


MOlO 


-19.5 


Lll 


8.41 


HERA 


-0.01 


HERA 


NGC 3368 


10.52 


Lll 


9.00 


MAIO 


-20.0 


Lll 


8.30 


Y95 


-0.45 


CIO 


NGC 3521 


10.70 


W08 


8.70 


MOlO 


-20.3 


Lll 


8.96 


HERA 


0.34 


HERA 


NGC 3627 


9.30 


W08 


8.66 


MOlO 


-20.1 


Lll 


8.84 


HERA 


0.36 


HERA 


NGC 3628 


9.40 


Lll 


9.00 


MAIO 


-19.6 


Lll 


9.20 


Y95 


0.33 


CIO 


NGC 3938 


12.20 


W08 


8.70 


E08; CIO 


-19.6 


W08 


8.41 


HERA 


-0.07 


HERA 


NGC 4194 


42.00 


CIO 


8.70 


E08; CIO 


-20.5 


NED 


8.90 


Y95 


1.13 


CIO 


NGC 4254 


20.00 


W08 


8.79 


MOlO 


-21.3 


W08 


9.50 


HERA 


0.83 


HERA 


NGC 4321 


14.30 


W08 


8.84 


MOlO 


-20.9 


W08 


9.20 


HERA 


0.45 


HERA 


NGC 4449 


4.20 


W08 


8.30 


M97 


-18.1 


Lll 


7.01 


B03 


-0.45 


CIO 


NGC 4450 


27.10 


CIO 


8.90 


CIO; MAIO 


-21.7 


NED 


8.90 


Y95 


-0.18 


CIO 


NGC 4536 


14.50 


W08 


8.60 


MOlO 


-19.7 


W08 


8.60 


HERA 


0.42 


HERA 


NGC 4569 


20.00 


W08 


8.90 


E08; CIO 


-22.1 


W08 


9.14 


HERA 


0.29 


HERA 


NGC 4579 


20.60 


W08 


9.00 


CIO; MAIO 


-21.4 


W08 


8.94 


HERA 


0.11 


HERA 


NGC 4631 


8.90 


W08 


8.43 


MOlO 


-19.9 


Lll 


8.72 


HERA 


0.40 


CIO 


NGC 4725 


9.30 


W08 


8.73 


MOlO 


-20.2 


W08 


7.85 


HERA 


-0.43 


HERA 


NGC 4736 


4.70 


W08 


8.66 


MOlO 


-19.4 


Lll 


8.14 


HERA 


-0.29 


HERA 


NGC 4826 


7.50 


W08 


8.87 


MOlO 


-20.0 


Lll 


8.10 


H03 


-0.50 


CIO 


NGC 5033 


14.80 


MOlO 


8.66 


MOlO 


-20.8 


NED 


9.30 


H03 


0.10 


K03 


NGC 5055 


10.10 


W08 


8.77 


MOlO 


-20.7 


Lll 


9.10 


HERA 


0.34 


HERA 


NGC 5194 


8.00 


W08 


8.86 


MOlO 


-20.6 


Lll 


9.20 


HERA 


0.49 


HERA 


NGC 5236 


4.50 


W08 


9.00 


MOlO 


-20.1 


Lll 


8.90 


Y95 


0.37 


CIO 


NGC 5253 


3.15 


Lll 


8.20 


MAIO 


-16.6 


Lll 


5.80 


T98 


-0.22 


CIO 


NGC 5713 


26.50 


W08 


8.64 


MOlO 


-20.9 


W08 


9.17 


HERA 


0.76 


HERA 


NGC 5866 


15.10 


CIO 


8.70 


CIO; MAIO 


-20.2 


NED 


8.10 


Y95 


-0.60 


CIO 


NGC 5953 


35.00 


CIO 


8.70 


E08; CIO 


-20.0 


NED 


9.00 


Y95 


0.38 


CIO 


NGC 6822 


0.49 


GIO 


8.11 


LE06 


-15.2 


NED 


5.15 


GIO 


-1.85 


Ell 


NGC 6946 


5.90 


W08 


8.73 


MOlO 


-19.2 


Lll 


9.04 


HERA 


0.57 


HERA 


NGC 7331 


14.70 


W08 


8.68 


MOlO 


-21.7 


NED 


9.10 


HERA 


0.49 


HERA 



^ References: B03 = B ottner et al.l ((20031 ): CIO = ICalzetti eTall f2010); D84 = IDufourl (fT98l ): EOS = lEngelbracht et al.l 
(2008); Ell = Efremova et al.l (120111); F08 =^ Fukui et al. (2008); G03: Gil de Paz et al. (2003); G04: Galleti et al. (2004); G08 
= Gro cholski et al. (2008^): GIO = IGratier et al. (2010a); HOI = Hunter (2001); H03 = Heifer et al. (2003); H04 = Hever et aQ 
(|200^ ) - H09 = Harris & Zaritskv (2009); 197 = Israel (1997); 104 = Izotov & Thuan (2004); K03 = Kennicutt et al. (200|)" 
L79 ^iLeoueux et al. (1979); L03 = Lee et al. (2003); LE06 = Lee et al. (2006); L06 = Lerov et al. (2006); L07 = Lerov et al] 
{2007"); Lll = Lee et al. (2011); M97 = Martin (1997); M06 = Mizuno et al. (2006); MAIO = Marble et al. (2010); MOlO = 
[Moustakas et al. (2010); Pll = Pasquali et al. (2011); R08 = Rosolowskv & Simon (2008); T98 = Taylor et al. (1998); T05 
= iThuan fc Izo tov (2005); W04 = iWilke et al.i (i2004 ): W08 = iiWalter et al.i (i2008i ): Y95 = lYoung et al.i (il995 V: HERA = 
HERACLES collaboration. 
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Table 3 




CO Luminosities of Molecular Clouds 


Name 


^CO 1-0 
(K km s-i pc2) 


Reference 


M 33 EPRB 1 


1.8 X 10^ 


Rosolowskv et al. f2003) 


LMC N197 


7.0 X 10^ 


Fukui et al.,, (2008) 


SMC N84 


1.3 X 10^ 


Mizuno et al. (2001) 


ICIO Blla 


7.6 X lO'^ 


Lerov et al. (2006) 


Orion-Monoceros 


8.6 X 10^ 


Wilson et al. (2005) 


Orion A 


2.7 X lO'^ 


Wilson et al. (2005) 


Taurus 


5.6 X 103 


Goldsmith et al. (2008) 



(i.e., map coverage), and stacked over regions bright in 
24/im. 

3.1. Individual Lines of Sight 

We start with searching for significant CO emission in 
individual hnes-of-sight. For the dwarf galaxies in HER- 
ACLES the noise per channel map in the full resolution 
(13'' X 2.6 km s"^) cubes is a = 21 ± 3 mK. For each 
galaxy we search the entire cube for regions with signal- 
to-noise ratio (SNR) > 4 over two consecutive velocity 
channels. This corresponds to a CO point source with 
luminosity Leo 2-1 = 2.1 x 10"^ (cr/20 mK) {D/4 Mpc)^ 
K km s~^ pc^. For comparison, Table[3]lists CO(l-O) lu- 
minosities of the brightest clouds in M 33, LMC, SMC, 
IC 10, and values for the Milky Way Orion-Monoceros 
complex, Orion A, and Taurus. We are sensitive enough 
to detect these clouds (except Taurus) at our average 
source distance of D = 4 Mpc. 

Figure [2] shows maps of integrated CO intensity for the 
more massive dwarf galaxies of our sample. Each line-of- 
sight integral includes all channels with velocities within 
the range of ±50 km s~^ of the local Hi mean velocity. 
Five galaxies, NGC 2403, NGC 2976, NGC 3077, NGC 
4214, and NGC 4625, show emission exceeding our point 
source sensitivity within ±50 km s~^ of the local mean 
Hi velocity. A point source of 1.5 times our point source 
sensitivity will show up completely black at the chosen 
linear grayscale. For all other galaxies we detect no signal 
at this angular resolution. The non-detection of bright 
CO clouds in most of our targets is most likely linked 
to their low metallicity, 12 ± logO/H < 8.0, while the 
reference sample in Table [3] has higher metallicities, 8.2 < 
12 + logO/H<8.8. 

3.2. Improve Sensitivity by Stacking 

The large map size of the HERACLES maps and the 
fine (13'') resolution as compared to the angular extent of 
the galaxies allow us to search for CO emission at many 
different locations inside the galaxies. We saw above that 
only a few galaxies have signal strong enough to be de- 
tected in individual lines-of-sight. Therefore, we now ap- 
ply the s tacking techni q ue de veloped and described in 
detail in ISchruba et al.l (|2011[ ). This method accounts 
for the velocity shift in the observed CO spectrum due 
to galaxy rotation or other bulk motion. This is done 
by re-adjusting the velocity axis of the CO spectrum of 
each line-of-sight such that the local Hi mean velocity 
appears at a common (zero) velocity in the shifted spec- 
trum. Under the assumption that the mean velocities of 
Hi and CO closely correspond to each other (confirmed 
in the bright inner disk of spiral galaxies), the CO line 



peaks in the shifted spectra by construction at zero ve- 
locity across each galaxy (and across the sample). By 
averaging these shifted spectra we can decrease the noise 
and coherently add up the spectral line at known (zero) 
velocity. We may expect that the signal-to- noise (SNR) 
in the integra ted line intensity, Leo, does improve pro- 
portional to v^AT^^/v^Al^bs, where AVcai is the total 
velocity gradient across the galaxy from galaxy rotation 
and AVbbs is the width of the CO line at the observ- 
ing resolution. This ratio may typically be on the order 
of V250/a/20 ~ 3.5. Significant further improvement in 
the SNR will be achieved by averaging over many lines- 
of-sight. 

To determine the line intensity we fit the stacked spec- 
trum by a Gaussian profile with center restricted to be 
within ±50 km s~^ of zero velocity, FWHM to be larger 
than 15 km s~^, and the amplitude to be positive. In 
cases where the fitted Gaussian has peak intensity be- 
low 3<j or the integrated intensity is less than 5 times its 
uncertainty we determine an upper limit instead. The 
upper limit is defined as the integrated intensity of a 
Gaussian profile with FWHM set to 18 km s~^ and am- 
plitude fixed to 3a. Figure [3] shows stacked CO spectra 
determined over the entire galaxy for the targets shown 
in Figure [2l 

For white noise <Jrms (in the velocity-integrated inten- 
sity of the stacked spectrum) decreases proportional to 
A^~^/^ by stacking where TV is the number of indepen- 
dent resolution elements. For our data <Jrms does improve 
by stacking but at somewhat slower rate and saturates at 
~ 1 mK km s~^ after averaging over ~ 500 — 1000 resolu- 
tion elements. Deviations from the white noise behavior 
are linked to our obse rving strategy and data reduction 
fsee lLerovet al.ll2009f ). 

It is further instructive to note that the size of the 
selected stacking region becomes a critical quantity as 
we have to deal with non-detections. Increasing the size 
of the stacking region will lead to stronger upper limits 
on /co; in an ideal world Iqo oc N~^^'^. The upper limits 
on Leo will however degrade with increasing area since 
Leo = ^co X Area ex A^^/^. As we are mainly interested 
in the absolute quantity Leo, we have to be careful when 
selecting an appropriate stacking region. 

3.2.1. Stacking of Entire Galaxies 

We start with stacking the CO spectra over the entire 
(mapped) extent of each galaxy. To define the "entire" 
galaxy extent, we use the SFR distribution as guideline 
for the (most likely) distribution of molecular gas and CO 
emission. Unfortunately, this method does not provide 
definite sizes. FUV emission (the main tracer of SFR in 
dwarf galaxies outside massive star- forming regions) typ- 
ically starts to fiatten as function of galactocentric radius 
before reaching the background level. We therefore select 
for each galaxy a maximum galactocentric radius (typi- 
cally between 1 — 2 R25) that includes most (~ 95%) of 
the galaxy-integrated star formation. The so selected re- 
gions are highlighted by black contours in Figures [T] & [2l 

We will show later that for some galaxies the as above 
selected region is not fully sampled by our CO map 
and may miss a significant fraction (in the worst cases 
- 10 - 30%) of the total SFR as given in Table [H This 
is especially true for NGC 5474 and NGC 2403, and to 
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Figure 3. Stacked CO Spectra. The resulting mean CO(2-l) spectrum after stacking all data over the entire galaxy (or map extent; black 
contours in Figures [T] &; (2) • The spectra are shifted to the local mean Hi velocity and thus expected to peak at zero velocity (see text). 
We fit Gaussian profiles (green lines) to the data within ±50 km s"-*^ to determine the integrated CO intensity; for stacked spectra without 
robust signal we determine Sa upper limit (green dot-dashed lines). The horizontal dotted lines show the la rms noise of the stacked 
spectra. 
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Table 4 

Stacking of Entire Galaxy 



Table 5 

Stacking of 24/im-bright Regions^ 



Name 



Area^ 

(kpc2) 



(10^ K km s-i pc2) 



SFR^ 

(M0 yr-i) 



Name Area i^24/2m^ Leo 2-1 SFR 

(kpc^) (10^ MJy sr-1 pc^) (10^ K km s'^ pc^) (Mq yr"!) 



M81 DwA 


1.33 


< 0.38 


0.0004 


M81 DwA 










M81 DwB 


1.79 


< 0.64 


0.0022 


M81 DwB 










DDO 053 


1.75 


< 0.37 


0.0035 


DDO 053 


0.31 


0.20 


< 0.16 


0.0016 


DDO 154 


7.74 


< 1.00 


0.0054 


DDO 154 










DDO 165 


11.47 


< 1.46 


0.0093 


DDO 165 










HO I 


12.37 


< 1.81 


0.0093 


HO I 










HO II 


26.69 


< 2.83 


0.0379 


HO n 


1.82 


1.35 


< 0.38 


0.0136 


IC 2574 


83.85 


< 8.20 


0.0670 


IC 2574 


2.39 


1.54 


< 0.52 


0.0158 


NGC 2366 


12.23 


< 0.86 


0.0532 


NGC 2366 


2.80 


7.25 


< 0.39 


0.0403 


NGC 2403 


39.40 


26.79 ± 0.31 


0.3131 


NGC 2403 


33.01 


50.87 


26.04 ± 0.28 


0.3021 


NGC 2976 


19.21 


13.89 ± 0.29 


0.0908 


NGC 2976 


7.63 


17.78 


12.10 ± 0.16 


0.0856 


NGC 3077 


7.78 


3.54 ± 0.15 


0.0689 


NGC 3077 


6.13 


19.10 


3.53 ± 0.13 


0.0679 


NGC 4214 


14.75 


3.21 ± 0.25 


0.1041 


NGC 4214 


5.99 


15.73 


2.33 ± 0.12 


0.0934 


NGC 4236 


59.54 


< 3.53 


0.1063 


NGC 4236 


7.15 


5.75 


< 1.01 


0.0500 


NGC 4625 


31.19 


5.73 ± 0.65 


0.0664 


NGC 4625 


9.48 


10.33 


5.12 ± 0.31 


0.0592 


NGC 5474 


39.19 


< 4.92 


0.0683 


NGC 5474 


7.64 


3.16 


< 2.19 


0.0375 



^ Unprojected area sampled from this work. 
^ -f/CO(2— 1) i^ sampled area from this work. 
^ SFR(FUV+24) in sampled area from this work. 

a lesser extent for Ho II, NGC 2366, NGC 3077, NGC 
4214, and NGC 4236. It is, however, not obvious how 
to correct for this effect. In the remainder of this paper 
we will therefore continue to refer to our measured Leo 
as the total galaxy-integrated luminosity but urge the 
reader to keep in mind that the true value may be up to 
~ 30% higher for a (small) subset of our sample. The 
given uncertainties on Leo include only the statistical 
uncertainties of fitting the stacked spectrum with Gaus- 
sian profiles. Uncertainties in the calibration (from in- 
strumental and re duction met hodology) may affect Leo 
by up to 30% (see lLerov et al.l 2009) and uncertainties in 
the distance will enter quadratically — neither effect is 
included. 

Figure [3] shows the resulting spectra when stacking over 
the entire observed part of the galaxy (for the same galax- 
ies shown in Figure [2j). Table |4] lists the (unprojected) 
area and the respective Leo and SFR. Five galaxies, 
NGC 2403, NGC 2976, NGC 3077, NGC 4214, and NGC 
4625 are robustly detected. These are the same galax- 
ies that already showed emission for individual lines-of- 
sight (Section [5TT]) . One galaxy, NGC 4236, may show a 
tentative signal which extends from —8 to +20 km s~^, 
has peak intensity ~ 2.4 mK (~ 2.3<j) over 2 channels, 
Ico 2-1 ~ 0.035 K km s-\ and Leo, 2-1 ~ 2.1 x 10^ 
K km s~^ pc^, a factor 0.6 below our quoted upper limit. 
This emission is not point-source-like because with a 
point source sensitivity of Leo ~ 2.5 x 10^ K km s~^ pc^ 
for this galaxy it would have been easily detected. All 
other galaxies remain undetected. 

There are three galaxies, IC 2574, Ho II, and NGC 
5474, where we may have expected to find signal as 
these galaxies have properties similar to detected galax- 
ies. The stacked spectrum of IC 2574 shows an enhance- 
ment peaking at ~ — 5 km s~^, with full width ~ 16 — 18 
km s~^, and peak intensity ~ 4 mK (~ la) for the 24/im- 
selected regions or ~ 1.3 mK (~ 0.8<j) over the entire 
galaxy. While the match between CO and Hi velocities 
is encouraging, the significance of this enhancement is 
too low to diflFerentiate it from spurious emission. The 
stacked spectra of Ho II and NGC 5474 show no signs of 



^ This region includes all lines-of-sight with /24/xm > 0.2 MJy sr 
at 13^^ resolution. 

^ These units allo w comparison to the trend /24/xm ~ Ico found 
for massive spirals (jSchruba et al.ll2011l ). 

signal at a noise level of 1.9 and 2.2 mK per 2.6 km s~^ 
channel, respectively. 

3.2.2. Stacking of 24 l^m-bright Regions 

We make a final attempt to search for faint CO emis- 
sion by stacking over regions that likely have the highest 
probability to contain molecular gas and may be bright 
in CO. These are regions rich in dust and showing signs 
of embedded high-mass star formation. We use the 24/im 
intensity, /24/im, as a tracer of these conditions and select 
all lines-of-sight that have /24/im ^ 0.2 MJy sr~^ at 13^' 
resolution. The adopted 24/im level does not have a spe- 
cific physical interpretation, but it is well (^ 4(j) above 
the noise level of the 24/im maps. IR emission tend s 
to be faint in dwarf galaxies (e.g. IWalter et al.l 120071 ). 
In the more massive dwarf galaxies of our sample this 
24/im level effectively separates the star-forming peaks 
from the rest of the galaxy. The smallest dwarf galaxies 
however do not reach this 24/im intensity, and we omit 
them from this analysis. The thus selected regions are 
highlighted by red contours in Figure [TJ 

Table miists the results when stacking over these 24/im- 
bright regions: the (unprojected) area, 24/im and CO lu- 
minosity, and enclosed SFR. This method does not lead 
to new CO detections in addition to those galaxies al- 
ready detected at individual lines-of-sight and over the 
entire galaxy. However, for non-detected galaxies it re- 
sults in stronger upper limits on integrated quantities 
that scale with the size of the stacking regions, i.e., lower 
upper limits on Leo and lower Leo /SFR ratios. 

4. SCALING RELATIONS FOR CO LUMINOSITY 

4.1. Comparison to Magellanic Clouds 

We begin with a comparison of our CO measurements 
(upper limits; Table |4] & [5j) to the Magellanic Clouds 
(Table [2j). These are essentially the only low-metallicity 
systems that are well detected in CO over the full galaxy 
extent. We also list previous CO observations of our tar- 
gets in Table [6l A direct comparison to our CO measure- 
ments is however not straightforward as previous obser- 
vations covered only small fractions of the star-forming 
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Table 6 

Previous CO Observations of our Galaxy Sample 



Name Beam ^co i-o^ 

(arcsec) (10^ Kkms"-*^ pc^) 



Reference 



M81 DwA 


45 


< 0.16 


M81 DwB 






DDO 053 


55 


< 0.94 


DDO 154 


65 


< 3.2 


DDO 165 


65 


< 5.6 


Hoi 


65 


< 3.8 


Ho II 


13x65 


< 9.9 




45 


< 0.24 




55 


< 1.2 


IC 2574 


10x65 


< 10 




55 


1.1 


NGC 2366 


3x65 


< 1.9 




22 


< 0.50 




22 


< 0.04 




55 


< 1.3 


NGC 4214 


45 


0.38 




60 


0.56 




4x55 


0.73 


NGC 4236 


11x45 


< 9.4 


NGC 4625 


22 


4.4 




22 


4.2 




55 


16 


NGC 5474 


55 


< 1.5 



Young et al. (1995) 

Ler ov et al. (2005) 
Mor ris fc Lo (1978) 
Tavlor et al. (1998) 
Tavlor et al. (1998) 
Elmegreen et al. (1980) 
Young et al. (1995) 
Lerov et al. (2005) 
Elmegreen et al. (1 980) 
Lerov et al. (2005) 
Elmegreen et al. (1980) 
Hunter &: Sage (1993) 
Albrecht et al. (2004^ 



Lerov et ^l^ 
Young et al. 
Israel (19971 
Tavlor et al. (1998) 
Young et al. (1995) 
Boker et al. (2003) 
Albrecht et al. (2004) 
Lerov et al. (2005) 
Lerov et al. (2005) 



^ Luminosities are given on main beam temperature scale (Tmb) 
and are calculated assuming our adopted distances; line width and 
upper limit (3-4cr) on peak intensity are taken from the relevant 
reference. 

disk (often just a single pointing) and were strongly lim- 
ited by sensitivity. The large scatter between individual 
literature measurements and compared to our values in- 
dicate that previous pointed observations have not been 
able to robustly constrain the galaxy-integrated CO lu- 
minosity of dwarf galaxies. 

The galaxies that we detect in CO are comparable to 
(or exceed) the LMC in Mhi, Mb, SFR, and metallicity, 
but we are able to detect them at distances D = 2.9 — 9.5 
Mpc. Galaxies that have not been detected when stack- 
ing over the entire galaxy extent, have Leo 2-1 upper 
limits 0.1 — 2.6 times the CO lumi nosity of the LMC, 
^co'i-o = 3.2 X 10^ K km s"^ pc^ (iFukui et al.| [2QQ8). 
Our data is not sensitive enough to detect a CO lumi- 
nosity comparable to the SMC, ^^^^^^-0 = l-^ x 10^ 
K km s~^ pc^ (jMizuno et al.l [2QQ6f ) if it is spread over 
many resolution elements. For the IR-selected regions, 
our CO sensitivity improved and is always sufficient 
to detect Lqq^^_q and reaches down to 1 — 10 times 
^co^i-o- -^^^ individual lines-of-sight we would have eas- 



ily detected L^q^i-o t>ut detect no such point sources for 
11 of our 16 galaxies. 

4.2. Scaling Relations of Lqo with Mb and Metallicity 

We use our robust estimates of the galaxy-integrated 
CO luminosity of dwarf galaxies to examine the rela- 
tionship between Leo 1-0, 5-band magnitude, Mb, and 
metallicity in Figure ID In conjunction with our liter- 
ature compilation our galaxy sample covers 5 orders of 
magnitude in Leo, 4 orders of magnitude in Lb, 1.5 or- 
ders of magnitude in metallicity, 5 orders of magnitude in 
star formation rate (SFR = lO""^ - 10^ M© yr"^), and 3 
orders of magnitude in H I mass (Mhi = 10^ — 10^^ Mq). 

The top panel of Figure |4] shows Leo as a function of 
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Figure 4. Scaling Relations for CO Luminosity. Galaxy- 
integrated CO (1-0) luminosity as function of B-band magnitude 
(top) and metallicity (middle), and the ratio of Lqo and Lb as 
function of metallicity (bottom). Bigger symbols show stacking 
results of this work, smaller symbols show a compilation of liter- 
ature measurements. Color highlights metallicity. The solid line 
in the top panel shows a constant Lco/L-q ratio intersecting the 
bright galaxies. The trends are correlated by the well established 
luminosity-metallicity relation. 
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Mb- For guidance we show the sohd hne which high- 
hghts a constant scahng between Leo and Lb, set to in- 
tersect the bright galaxies. In the bright galaxies (Mb < 

— 18), Leo and Mb track one another with a more- 
or-less fixed ratio, logio Lco[K km s~^ pc^] /Lb[L0] = 

— 1.5 ±0.4, as one mi^ht expect for a s imple scaling with 
galaxy mass (I Young fc Scovilld [TOOTl iLerov et al.l I2QQ5I : 
iLisenfeld et"aLll2Qll[ ). The dwarf galaxies (Mb > -18) 
on the other hand lie below the solid line. They are 
"underluminous" in CO, i.e., their ratios are systemati- 
cally smaller, — 2.7 ± 0.6, than those of massive galaxies. 
Despite this trend, Leo and Mb are strongly correlated 
with (absolute) rank correlation coefficient rcorr = 0.79. 

The middle panel of Figure H] shows Leo as function 
of metallicity. There is a dramatic drop in Leo by 3 — 4 
orders of magnitude over a small range of metallicities. 
This drop is to first order caused by the much smaller 
mass and size of dwarf galaxies. However, due to the 
strong lu minosity-metallicity relation for dw arf irregu- 
lars fe.g.. iLee etaH 120061 : iGuseva et al]l2009[ ). it is also 
correlated to Mb (i.e., the top panel). The rank correla- 
tion coefficient is rcorr = 0.60. 

The bottom panel of Figure |4] shows the ratio Lqo/Lb 
as function of metallicity. Plotting Leo/ Lb should re- 
move most of the mass and size dependence seen in the 
above panels. The decreasing trend of Lqo/Lb with de- 
ceasing metallicity clearly shows that dwarf galaxies are 
also "underluminous" in CO in a normalized sense. 

From the s tudy of a large sample of literature CO data, 
iTavlor et al.l (|1998 ) suggested a "detection threshold" for 
CO below 12+logO/H ^ 8.0, about the metallicity of the 
SMC. Our data do not overcome this threshold: All of 
our galaxies with lower metallicity remain undetected. 
However, given the decreasing trend of Leo /Lb with de- 
creasing metallicity and the fact that our data is not 
sensitive enough to detect a galaxy like the SMC at a 
distance of D = 4 Mpc leaves open the question whether 
the proposed threshold is of observational or physical ori- 



gin. 



4.3. CO and Tracers of Star Formation 



Figure [5] shows the relationship between CO emission 
and tracers of recent star formation, FUV and 24/im, to- 
gether with a combination of FUV and 24/im often used 
to estimate the SFR (see Section [23|) . Panels in the left 
column show the correlations between observables /eo, 
/fuv7 and /24/im , panels in the right column show the ra- 
tios of /eo and /fuv or /24/xm as function of metallicity. 
For our dwarf galaxies (bigger symbols) we show both 
the stacking results derived over the entire galaxy (col- 
ored by metallicity) as well as the values derived from 
stacking over the 24/im-bright regions (gray symbols). 
For comparison we also show azimuthally averaged radial 
intensity profiles for spiral g alaxies (smaller colore d sym- 
bols) which are taken from iSchruba et ahl ()2011f ). The 
plots shown here are similar to plots of the "star forma- 
tion law", i.e., plots of Ssfr versus I;h2, though these 
are typically presented with axes interchanged and show 
data that are corrected for inclination (which we have 
not done here). 

One of the results from ISchruba et al.l (|2011[ ) was that 
spiral galaxies exhibit a strong correlation between /eo 
and /24/im (see middle panel in left column). This has 
been interpreted as a direct link between molecular gas as 



traced by CO emission and SFR, which is mostly deeply 
embedded and traced by 24/im emission. The ratio of 
CO/24/im is roughly constant inside galaxies and shows 
only little variation between galaxies and as function of 
metallicity where CO is detected (middle panels). The 
ratio of CO/FUV also shows little scatter inside individ- 
ual galaxies but can vary significantly between galaxies 
(upper panels). The large scatter in CO/FUV and the 
strong scaling with metallicity refiects the strong sus- 
ceptibility of FUV emission to dust attenuation. In low 
metallicity environments CO/FUV is low because CO 
abundance and thus /eo is low and at the same time, 
low dust abundance and low attenuation cause /fuv to 
be relatively high. The ratio of /eo/(/FUV+0.04/24/xm) 
which is proportional to the H2 depletion time, Tdep, 
is to first order constant ~ 1.8 Gyr for environments 
with 12+logio 0/H ^ 8.7 and shows little dependence on 
metallici ty in environments with about solar meta llicity 
(see also lBidel et~aII[2QQ8ll2QTTI : ILerov et al.ll201lD . 

Our dwarf galaxies however do show deviations from 
the trends observed for radial profiles of more massive 
galaxies. The CO, FUV, and 24/im intensities measured 
in the dwarf galaxies are close to the lowest intensities 
measured in the radial annuli of more massive galaxies. 
In addition, the ratios of CO/FUV and CO/24/im are 
shifted to smaller values. For our detected galaxies the 
ratios are a factor 5 — 10 below the ratios found in more 
massive galaxies. The data of our undetected galaxies are 
scattered but for galaxies with sensitive CO upper limits 
they are also shifted toward low CO/FUV and CO/24/im 
ratios. Dwarf galaxies exhibit enhanced signatures of star 
formation (both embedded and unobscured) per unit CO 
brightness as compared to large star- forming galaxies. 

5. IMPLICATIONS FOR CO-TO-H2 CONVERSION FACTOR 

A serious complication in studying the molecular con- 
tent of dwarf galaxies arises in how to relate the ob- 
served CO luminosities to H2 masses. Applying a 
Galactic CO-to-H2 conversion factor, aeo,Gab to dwarf 
galaxies that h ave been d etected in CO results in 
low H2 masses jTavlor et al.l 119981 : Mizuno et alj 120011 : 
iLeroy et al._ 20071 ). The resulting H2 masses are so low 
that to explain the observed SFRs the conversion of H2 to 
stars would need to be on average 10 — 100 times more ef- 
ficient than in Galac tic environments — a condition that 
seems unlikely (e.g.. iBolatto et ani2011f ). 

The detection of excess ionized carbon and infrared 
to millimeter dust emiss i on around star- forming; re- 



gions (iMadden et al .iri997l : iPak et al.lll998l: iRubin et al.l 
120091 : iCormier et al.l 120101 : llsrael fc Malonevll2011f ) indi- 



cate s that CO may not tra ce all H2 at low metalhc- 
ity ([ A MoiLeVfc Black 19881 llsraell 119971 : [Bolatto et^ 
119991 : IWolfiree t al. 2010). Because H2 can self-shield, 
its abundance is basically a function of its formation 
time (which depends on metallicity), however, CO can- 
not self-shield and exists only in regions that are suffi- 
ciently shielded by dust from the interstellar rad iation 
field (Glover et al. '20T0!: IGlover fc Mac L^[20TTh . o^eo 
is therefore assumed to be a strong function of metallicity 
and radiation field strength, although robust functional 
paramet rizations of these d epen dences are still lacking 
(but see lShettv et al.ir2011al lbrand lNaravanan et al.ll201ll 
for recent theoretical works). 
In the following we will discuss three different meth- 
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Figure 5. CO Emission and Tracers of Star Formation. Left column: CO (1-0) intensity as function of FUV and 24/im intensity, and 
a combination of FUV and 24/i,m used to determine the SFR. Right column: Intensity ratios of CO and FUV or 24/im as function of 
metallicity. Bigger symbols show stacking results for dwarf galaxies derived over the entire galaxy extent (colored symbols) or over 24yLim- 
bright regions (gray symbols). Smaller symbols show azimuthally averaged radial profiles in massive spiral galaxies. CO emission is well 
correlated with emission of SFR tracers, especially with 24/im. The ratios CO/FUV and CO/24/im show systematically smaller values in 
low-met allicity environments as is typical for dwarf galaxies. 
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ods that have been appHed to estimate aco in external 
galaxies. In particular we are interested in the metallicity 
dependence of aco- Whenever possible, we parametrize 
this dependence by a linear function between logio o^co 
and (12+logO/H), i.e., 

logio c^co = logio A + TV X (12 + logio 0/H - 8.7) (1) 

where the normalization A gives aco at 12+logO/H = 
8.7 and N is the slope. 

Before we begin, we have to caution the reader that gas 
phase metallicities bear considerable uncertainties. Dif- 
ferent empirical and theoretical calibrations can result 
in systematic discrepa ncies in estimate d metallicities as 
large as 0.1-0.7 dex (Kewley & Ellison 2008). However, 
once a specific calibration is selected the relative order- 
ing of individual galaxies and derived slopes are more 
robust. For our new data and our literature compilation 
we have tried to maximize homogeneity of metallicity es- 
timates (i.e., we use the metallicity calibration described 
in Section 12.41 whenever possible). To be on the save 
side, we thus concentrate our discussion on the slope of 
the parametrization between aco and metallicity. For 
the same reason a direct comparison of our results to the 
literature is hindered as previous studies applied a va- 
riety of different metallicity calibrations that may even 
vary within individual studies. 

5.1. Different Methods to Estimate the CO-to-H2 
Conversion Factor 

5.1.1. Virial Method 

The classic method to derive aco uses high resolution 
CO observations that are capable of resolving individ- 
ual molecular clouds (e.g., [Solomon et al.l 119871 ). Under 
the assumption that a CO-bright core is in virial equi- 
librium, its observed linewidth and size can be converted 
into a virial m ass, Myi r . and from that aco = Myir/ Lco- 
Early work bylWilsonl (|T995h . I Arimoto et al.l (|T996l ). and 
iBoselli et al.l (|2002[ ) have applied this method to a hand- 
ful of Local Group galaxies and found a weak metal- 
licity dependence of aco with slopes flatter than —1. 
This metallicity dependence however has no t been con- 
firmed by the recen t studies of .Blitz et al.l (|2007[ ) and 
iBolatto et al.l (|2008[ ). They re-analyze a large set of 
literature data aiming at maximizing homogeneity of 
their analysis and carefully correcting for finite spa- 
tial and spectral resolution. They derive a distribution 
of aco values that scatters without systematic trend 
around 0.5 — 5 aco,Gai (the green striped region in Fig- 
ure [71 we describe this figure later in detail). It has 
to be emphasized that the virial method likely leads 
to a significant underestimate of the total H2 mass of 
GMCs in low-metallicity environments. This bases on 
the idea that with decreasing metallicity the CO-bright 
cores shrink to ( much) smaller size than the surrounding 
H2 clouds (e R., IBolatto et al.l 119991 : lWolfireet~aII 120101 : 
IShettv et al"][20m]lbl ). 

5.1.2. Dust Modeling 

This method uses IR observations and dust model- 
ing to estimate the gas mass and distribution which 
has the advantag e that it is independent of CO emis- 
sion (fThronson eFan 119881 : [Israeli 119971 : iDame et al] 



[200TI: iLerov et al.ll 2007L [20Q1 [2QTTI : iGratier et al.ll2010bl : 
Bolatto et al. 2011). It builds on the assumption that gas 
and dust are well mixed and aco is derived from Mdust = 
DGR X (Mhi + aco^co), where DGR is the dust-to-gas 
ratio. By modeling the dust distribution the local H2 
mass can be inferred from Mdust (after subtraction of lo- 
cal Hi) by either fixing the DGR in quiescent, non-star- 
forming regions (assumed to be H2-free) or by simultane- 
ous optimizing aco and DGR such that the scatter be- 
tween Mdnst/D GR an d Mhi + aco^co gets minimized. 
Early work by [Israeli (^997) implied a strong metallic- 
ity dependen ce of aco with slope of —2.7 =b . 3. Re - 
cent work by IGratier etlill (|2010b[ ). iLerov et al.l (|2011[ ). 
and IBolatto et al.l ()2011[ ) did confirm a strong increase 
of aco at low metallicities, although their results vary 
in absolute terms, proposed functional form, and are 
systematically smaller than the aco values derived by 
Israel (1997). The aco values derived by these studies 
lie within the blue striped region in Figure [71 The lowest 
metallicity galaxy to which this method has been applied 
is the SMC, for which large amounts of H2 have been in- 
ferred implying aco values 10 — 100 times the Galactic 
value (Leroy et al. 2011; Bo latto et al. 2011). 

Disadvantages of this method are that it is suscepti- 
ble to variations of the FIR emissivity of dust grains and 
variations in DGR between dense, star-forming regions 
and low density, quiescent regions. There are indications 
that the em issivity is enh anced in dense re^ons (e.g . 
iParadis et al.l 120091 : IPlanck Co llaborat ion et all 120111 ) 
which would cause an overprediction of aco on scales 
of individual star-forming regions. Second, because dust 
enrichment of the ISM by stars seems insufficient to ex- 
plain observed dust abundances, it is proposed that most 
dust forms in the ISM, presumably in the densest re- 
gions (Dwek 1998; Draine et al. 20Q3). If this dust is 
only slowly transported into the lower density ISM then 
this would also lead to an overprediction of aco- The 
need for sensitive, matched high resolution data to make 
the analysis robust limits this method to nearby galaxies 
and makes observations time consuming. 

5.1.3. Constant SFE 

An alternative method to constrain the II2 mass is to 
assume that the conversion of H2 to stars is indepen- 
dent of environment, i.e., assuming a constant H2 deple- 
tion time, Tdep, or a constant star formation efficiency 
(SFE; the inverse of Tdep). <^co is then given by aco = 
'Tdep X SFR/Lco- This approach has been apphed even 
when it was still consider ed very uncertain how to relate 
CO to H2 in our Galaxy (|Rana fc WilkinsQn[ll986[ ). We 
will apply this method to our sample of nearby galaxies 
in the remainder of this paper. The idea is encouraged by 
several observations: (a) the accumulating evidence that 
star formation in molecular clouds is independent from 
environment as indicated by the similarity o f molecular 
cloud properties in our and nearby galaxies f Blitz et al.j 
120071 : IBolatto et al . 1 120081 : iFukui fc Kawamur a 2010), the 
universality of SFE per free-fal l time in clouds of differ- 
ent mass and density (Krumholz fc Tan 2007). and evi- 
dence in favor of a u niversal initial stellar mass function 
([Bastian et al.l [20101 ): and (b) the remarkably constant 
scaling between H2 and SFR on ~ kpc sca les observed 
in a large se t of nea rby spiral galaxi es (e.g., iBigiel et all 
[2OO8 .. aon,; ILerov et a h , 200& ; S chruba et al.l 1201 1[ ). A 
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Figure 6. Metallicity Dependence of the CO-to-H2 Conversion Factor, aco is derived from the ratio of observed SFR scaled by a 
constant H2 depletion time, r^ep = 1-8 Gyr, and the observed CO luminosity, Leo- Bigger symbols show galaxy-integrated measurements 
of dwarf galaxies from this work, smaller symbols show data for our literature compilation with starbursts highlighted by stars. The 
horizontal dashed line shows the Galactic conversion factor and the diagonal lines show regression fits: to all galaxies (dotted line), to all 
non-starburst galaxies (dashed line), and exclusively to the HERACLES sample (solid line). 



drawback of this method is that it requires that the cor- 
relation between H2 and SFR estabhshed in spiral galax- 
ies continues to hold in dwarf galaxies. This makes the 
method less rigorous than direct attempts to trace H2 
but also makes it available for a much larger sample of 
galaxies including distant galaxies. 

Currently we are not able to conclude that r^ep is 
truly constant. Observations of strongly variable star 
formation histories and starbursts readily indicate that 
it does no t hold in all environments (e.g. Lee et al."2009'; 
I Weisz et~a l . 2011). However , rece nt th eoretical consid- 
eratio ns by iKrumholz et al.l (|2Qllf ) and iGlover fc Clarkl 
(|2012[ ) provide a motivation in favor of a constant 
H2/SFR ratio. Although some of these methods do ques- 
tion if H2 is fundamental for star formation, they also ar- 
gue that H2 will be a good tracer of star-forming regions. 
This is because the Hi to H2 transition and the drop in 
gas temperature which makes clouds susceptible to grav- 
itational instabilities occur under similar conditions that 
are to first order set by dust shielding of the interstellar 
radiation field. 

5.2. New & Literature Measurements 

In the following we explore the implications for CO 
if Tdep is indeed constant. The average value for spi- 
ral galaxies in the HERACLES sample with about so- 
lar metallicity varies in the ran^e of ^ 1.8 to 2.35 Gyr 
([Bigiel et al. 2008 , 2QlTl lLerov et al.ll2QQ8l : lSchruba et al.l 
l2011f ). This natural variation is due to minor differences 
from galaxy to galaxy, and depends on exactly which 
lines-of-sight are included, and in particular the weight- 
ing method employed (i.e., galaxy average, radial rings, 
or pixel averages). We adopt here a value of Tdep = 1.8 



Gyr which is derived on the basis of entire galaxy aver- 
ages. Figure [6] shows the resulting aco values as func- 
tion of metallicity. In this plot we show galaxy-integrated 
values, bigger symbols show our measurements of dwarf 
galaxies (Table |4|) and smaller symbols show data from 
our literature compilation (Table [2]). Color coding high- 
lights metallicity as in previous plots. Star symbols in- 
dicate galaxies that are labeled in the literature as star- 
bursts. 

The derived aco values strongly depend on metal- 
licity. For galaxies with 12+logio 0/H > 8.6, we find 
aco ^ <^co,Gai although with ~ 0.3 dex (factor 2) 
scatter. For galaxies with lower metallicity, aco in- 
creases strongly with decreasing metallicity. For dwarf 
galaxies with 12+logio 0/H < 8.6, even though most 
of them remain undetected in CO, we can readily ex- 
clude aco ^ <^co,Gai- The few dwarf galaxies with CO 
detection suggest aco ^ 10 aco,Gai at 12+logioO/H 
< 8.4. We emphasize that the derived aco values for 
dwarf galaxies with Z/Zq ~ 1/2 — 1/10 are 1 — 2 orders 
of magnitude higher than aco values derived for massive 
spirals with Z/Zq ^ 1. 

We attempt to parametrize this dependence by fitting 
function of the form given in Eq. ([T]) . We use a bisecting 
linear regression to determine the best-fitting parame- 
ters. Uncertainties are determined from a Monte Carlo 
analysis. We have repeatedly added Gaussian noise to 
aco with log-normal standard deviation of 0.3 dex and 
to 12+logio 0/H with standard deviation of 0.1 dex and 
re-fitted the perturbed data. The quoted uncertainties 
correspond to the standard deviation of 100 such de- 
rived best-fit parameters. Table [71 lists the resulting 
normalizations and slopes together with the scatter of 
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Different Methods; 

CMC Modelling: 
Bolatto + 08 

Dust Modelltng: 
Leroy+ 1 1 

SFR Scaling: 

GcnzcM-li (z^l) 
This ^MoTk (z = 0) 




Table 7 

Metallicity Dependence of aco assuming a constant SFE^ 



Selected Data 



Value at 
12+logO/H=8.7 



Slope of 
Regression 



Scatter 
(dex) 



-2.8 ±0.2 
-2.4 ±0.3 



0.13 
0.10 



7.5 8.0 8.5 9.0 

Metalhcity [12 + Iogio0/H] 

Figure 7. Trends of the Metallicity Dependence of the CO-to-H2 conversion factor aco ■ Striped bands indicate the range of aco values 
derived from employing different methods (see text). The width in the bands indicate roughly the scatter of individual measurements. 

(dotted line). We consider the latter result uncertain and 
potentially biased high because it is driven by a handful 
galaxies that are currently undergoing a starburst and 
have CO measurements from Taylor et al. (1998), mea- 
surements that have been made prior to the latest gener- 
ation of sensitive millimeter receivers. The scatter of the 
data to the best-fit relations is 0.09 - 0.12 dex (- 30%) 
which is significant smaller than the scatter of '^ 0.3 dex 
in the ratio aco/o^co, Gai for galaxies with 12+logio 0/H 
> 8.6. In this sense a steep increase of aco with decreas- 
ing metallicity is much favored as compared to a constant 
value. The trend fitted to the "complete, non-starburst" 
sample and its associated uncertainty is indicated as red 

striped region in FigurejTj 

The recent study by Genzel et al.l (|2012l ) also applied 
the assumption of a constant SFE. They analyzed a 
sample of star-forming galaxies at redshift z ^ 1 — 2 and 
determined a slope of — 1.9 ± 0.67; the orange striped 
regiorrn in Figure [71 They also combined their distant 
galax y sample with the dust-inferred aco measurements 
from iLeroy etld] f201 1) which reduces their slope 
to — 1.3 =b 0.25. The decrease in slope is basically 
driven by two galaxies, M31 and the SMC, and their re- 
sult may be affected by combining two different methods. 



complete sample 

• all galaxies 8.2 ± 1.0 

• non-starbursts 6.9 ± 1.0 

HERACLES sample 

• all galaxies 8.0 ± 1.3 

• non-starbursts 7.1 ± 1.2 



-2.0 ±0.4 0.10 
-2.0 ±0.4 0.09 



^ Using bisecting linear regression of logio aco = logio A-\-N x 
(12±logio O/H — 8.7); see Eq. {T}. Uncertainties are determined 
by repeatedly adding random noise of 0.1 dex to 12 ± logio O/H 
(x-axis) and 0.3 dex to logio Q^CG (y-axis). 

the data orthogonal to the best-fit regression. We di- 
vide our galaxy sample in two groups: "starbursts" and 
"non-starbursts" , and evaluate the HERACLES galaxies 
and the complete sample separately. This may help to 
minimize biases due to inhomogeneous data sets. We ex- 
pect the smallest systematics for the HERACLES sample 
including only non-starburst galaxies. We separate the 
starbursts because they likely violate our assumption of 
a constant r^ep having SFR in excess of their H2 content. 
The best-fit regressions depend somewhat on the par- 
ticular galaxy sample; see Figures [6] and Table [3 For this 
analysis we neglect the upper limit measurements as they 
are not stringent enough to affect our best-fits. For the 
HERACLES sample we determine a slope of —2.0 ± 0.4 
roughly independent whether startbursts are included or 
not (solid line), but with larger uncertainties due to the 
relative small dynamic range sampled by the detected 
galaxies. For the complete galaxy sample, the slope is 
steeper. We determine a slope of —2.4 ± 0.2 for the non- 
starbursts (dashed line) and — 2.8 ± 0.2 for all galaxies 



5.3. Comparison 

Approximate trends for the metallicity dependence of 
aco derived from the three discussed methods and their 
intrinsic scatter are indicated in Figure [71 At solar metal- 
licities the three methods give roughly consistent results 



^'^ Note that we changed the normalizations of the Ge nzel et all 
parameterizations to match the data plotted in their Figure 3. For 
the high redshift sample we increased the normalization by 0.3 dex, 
for the combined sample we decreased it by 0.07 dex. 
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within their uncertainties. Toward lower metahicities the 
three methods however predict different trends for the 
dependence of aco on metahicity. The aco values de- 
rived from the virial method (green striped region) show 
no systematic trend with metallicity and exhibit roughly 
an order of magnitude scatter. The dust-inferred aco 
values (blue striped region) and the aco values derived 
from scaling the SFR (red and orange striped regions) 
however do show a strong systematic increase toward low 
metahicities. For dwarf galaxies with Z/Zq ~ 1/2 — 1/10 
the dust-inferred aco values are roughly a factor of 10 
larger than the Galactic value but with ~ 1 order mag- 
nitude scatter, for the SFR-scaled aco values they are a 
factor 10 — 100 larger. 

In addition to changes in the distribution of CO emit- 
ting gas at low metahicities, changes in the physical 
conditions may also affect the excitation temperature 
and thus the (relative) brightness of the CO (1-0) and 
CO (2-1) transitions. For molecular clouds in the Mag- 
ellanic Clouds, line ratios of R21 ^ 1.0 — 1.5 are fre- 
quent lyobserved (|Bolatto et al.ll2000L 120031 : llsrael et al.l 
I2OO3I . [2QQ5I ). This is typically (but not always) ascribed 
to higher fractions of warm and optically thin gas as CO 
is less shielded against the interstellar radiation field. In 
diffuse gas, values of R21 can be as high as ~ 3, how- 
ever, these diffuse regions are typically faint and con- 
tribu te only a minor f raction to the total CO emission 
(e.g., llsrael et"al]l2003f ). Thus our choice of using a sin- 
gle value of R21 = 0.7 when converting the HERACLES 
CO (2-1) line intensities to CO (1-0) values may bias our 
results to over-predict the true CO (1-0) line intensity in 
the HERACLES low metallicity galaxies by a factor < 2. 
The lack of a systematic offset to the sample of literature 
galaxies in Figures [6] that all use CO(l-O) data suggest 
that variations in R21 do not dominate our analysis. 

The discrepancy between aco values derived from the 
viral method and those derived from dust or SFR at low 
metahicities can be explained by considering the spatial 
scale which these methods operate on. The virial method 
operates on the small spatial scales of CO-bright cores 
of molecular clouds. At low metallicity these CO cores 
shrink as CO in low density gas gets dissociated while 
H2 can survive there via self-shielding. Applying such 
aco values to the total CO luminosity of a galaxy thus 
traces only the H2 mass within the high density CO- 
bright cores and at low metallicity will inevitably fail to 
trace the total H2 mass. On the other hand, the aco 
values derived from dust or SFR are sensitive to CO- 
dark H2 and can trace H2 on spatial scales of the size of 
star- forming regions and larger. 

But why are the aco values derived from dust and 
SFR different? We infer aco under the assumption of a 
fixed H2/SFR ratio, i.e., Tdep = 1.8 Gyr, that spans from 
massive spirals to low- mass, low- metallicity dwarfs. If 
this assumption breaks down then we will mis-attribute 
variations in H2/SFR to variations in aco- In models in 
which the star formation efficiency is set by the free-fall 
time, metallicity can affec t H2/SFR. T his is pointed out 
by Gnedin et aJ (2009) , Gnedin fc Kra vtsov (2011), and 
iFeldmann et al.l (|2011l ) ; they show that the gas densities 
containing H2 vary strongly with metallicity and radia- 
tion field, and thus free-fall times of clouds containing H2 
are not constant. While this can cause significant scat- 
ter in H2/SFR on cloud scales, on large (~ kpc) scales 



variations are expected to be much smaller. In environ- 
ments with metallicity Z/Zq ~ 1/10 and radiation field 
U/Uq - 10 - 100, free-fah times and thus H2/SFR are 
reduced by (only) a fa ctor 2 — 3. This is consistent with 
'Krumholz et al.l ( )2011l ): they find Ssfr/Sh2 to be con- 
stant within a factor 2 for ranges of I]h2 = 0.1 — 100 
Mo pc-2 and Z/Zq = 1 - 1/10. 

In the handful of studies that attempt to ac- 
count for aco variations an d measure Ho /SFR in 
low- metallicity galaxies (e.g., 'Gratier et ajj l2010al lbl: 
Bolatto et al. 2011), there are suggestions that H2/SFR 
is up to 2 — 5 times lower in local dwarfs at Z/Zq ~ 1/5. 
A factor of ~ 2 — 5 adjustment will not perfectly reconcile 
the various aco measurements at the lowest metahicities 
but can provide rough agreement at Z/Zq ~ 1/2 — 1/5. 
In this case Figure [71 and similar plots combine two im- 
portant trends: variations in aco and in Tdep- More 
detailed work comparing the SFR to H2 estimated via 
independent tracers like dust, [C II], or gamma rays will 
be needed to refine this approach. 

As noted above, the relative alignment of our galaxy 
sample along the x-axis (metallicity) in Figure [71 and pre- 
ceding plots are fairly secure. However the absolute cal- 
ibration of metahicities measured for extragalactic sys- 
tems remains uncertain. Therefore the relationship to 
solar metallicity remains somewhat tenuous, as does the 
alignment of the three methods. 

Another concern is that the ratios of CO, H2, and 
SFR vary with time. On small scales this shows 
as offsets between Ha, a tracer of recent star forma- 
tion, and CO, and induces large scatter in the re- 
spect ive ratios (e.g., iSchruba et al.ll2QlQl : lOnodera et al.l 
|2010[ ) which can be linked to the evol ution of star- 
forming regions (Kawamura e t all l2009f ). Numerical 
simulations even suggest that the ratios may never 
be constant inside a cloud because chemical equilib- 
rium is not reached during most or all of a cloud's 
lifetime (iGlover et al.l [20Tof iGlover fc Mac Lo^yl [Ml 

(|2009h 



Ti? 



'Shetty^eLa]] [2011al'b'; 'Feldmann fc Gnedin' "MAh.^ 
galaxy scales howev er Pelupessv fc Papadopoulos] ("" 
and Papadopoulos fc Pelupessy (2010) suggest that the 
ratios are roughly constant after dynamical equilibrium 
between ISM phases and stars is established {t ^ 1 Gyr). 
Only during (early) times of strong galaxy evolution, 
when the ISM phases and star formation are out of equi- 
librium, do larger deviations between CO, H2, and SFR 
occur {t < 0.2 — 0.3 Gyr); gas-rich and/or low- metallicity 
galaxies can show strong periodic variations through- 
out their evolution. Such variations are strongest in the 
smallest dwarf galaxies (Mb > —15) and are less com- 
mon and w eaker in (more) massive dwarfs and spirals 
(|Lee et al.l 2009). Some galaxies of our literature sample 
experience a current starburst (e.g., NGC 2366, NCG 
4449, and NGC 5253) or are in a post-starburst phase 
(e.^., NGC 1569). These bursts can last for a few 100 
Myr ( "McOuinn et al.|[2010[ ) and may show a strong time 
evolution in the brightness of their molecular gas and 
star formation tracers: starting with being bright in CO, 
followed by a phase being bright in CO and IR (a sign of 
embedded star formation), and finally being bright in IR 
and FUV (sensitive to stellar populations of age ^100 
Mvr: iSalim et al.]l2007[ ). Our literature samples, includ- 
ing SINGS, LVL, THINGS, and HERACLES have often 
an implicit or explicit bias to select IR-bright galaxies. 
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and thus actively star-forming systems which means that 
we infer a high aco- Robust volume-hmited surveys or 
otherwise unbiased samples are needed to remedy this. 

6. SUMMARY 

This paper presents sensitive maps of ^^CO J = 2 — 1 
emission for 16 nearby star- forming dwarf galaxies from 
the HERACLES survey (for a first presentatio n of a sub- 
sample of our galaxies see iLerov et alJ l2QQ9f ) . Thanks 
to the large area covered (^2 — 5 R25) and high linear 
resolution of ~ 250 pc at the average target distance of 
D = 4 Mpc, we can sample our targets by 10 — 1000 
resolution elements. 

We apply t he st acking techniques developed in 
iSchruba et al.l ()2011l ) to perform the most sensitive 
search for CO emission in low-metallicity galaxies across 
the entire star-forming disk. We search for CO emission 
on three spatial scales: individual lines-of-sight, stack- 
ing over IR-bright regions indicating embedded star- 
formation and thus regions likely to contain molecular 
gas, and stacking over entire galaxies. Our point source 
sensitivity is Leo 2-1 ~ 2 x 10^ K km s~^ pc^, suffi- 
cient to detect a CO-bright cloud with luminosity com- 
parable to Orion A or the brightest cloud in the SMC 
but at distance = 4 Mpc. When stacking over the 
entire galaxy our data have sufficient sensitivity to de- 
tect the LMC at i:> = 4 Mpc; but not the SMC. We 
detect 5 galaxies in CO with total CO luminosities of 
i^co2-i = 3 - 28 X 10^ K km s'^ pc^. The other 11 
galaxies remain undetected in CO even in the stacked 
images and have Leo 2-1 ^ 0.4 — 8 x 10^ K km s~^ pc^. 

We combine our dwarf galaxy sample with a large sam- 
ple of spiral galaxies from the literature to study the re- 
lations between Leo, ^b, and metallicity. We find that 
dwarf galaxies with metallicities Z/Zq ^ 1/2 — 1/10 have 
Leo of 2 — 4 orders of magnitude smaller than massive 
spiral galaxies with Z/Zq ~ 1 and that their Leo per 
unit Lb is 1 — 2 orders of magnitude smaller. Dwarf 
galaxies are thus significantly fainter in CO than a sim- 
ple linear scaling with galaxy mass would suggest. 

We also compare Leo with tracers of recent star for- 
mation (FUV and 24/im intensity) and find that Leo 
per unit SFR is 1 — 2 orders of magnitude smaller in 
dwarf galaxies as compared to massive spiral galax- 
ies. The low Leo /SFR ratios in dwarf galaxies may 
either indicate intrinsically small H2 masses coupled 
with high star formation efficiencies or that CO emis- 
sion becomes an increasingly poor tracer of H2. The 
two are degenerate, however, following the arguments of 
recen t observational studies of the dust-infe rred gas con- 
tent (Lerov e t al.]l2011l : iBolatto et all l2011') and theoret- 
ical studies o f the SFR-H2 dependence (Krumholz et al. 
IMlh Glover fc Clarkl l2012l we argue that the latter, i.e., 
significant changes in the CO-to-H2 conversion factor, 
Qfeo, at low metallicity are the dominant driver. 

To estimate aco and study its metallicit y dependence 
we ap ply a method recently also used by iGenzel et al.l 
(|2012f ) which assumes the conversion of H2 to stars to 
be constant and infer H2 masses and aco values by scal- 
ing the observed total SFRs. We assume an H2 deple- 
tion time of Tdep = Mh2/SFR =1.8 Gyr, the average 
valu e found for massive spirals in the HERACLES sam- 
pleJBigiel et al.l 120081 : ILerov et al] 120081 : ISchruba et al.l 
l2011f ). With this assumption we derive o^eo values for 



dwarf galaxies with Z/Zq ^ 1/2 — 1/10 more than one or- 
der of magnitude larger than those found in massive spi- 
ral galaxies with solar metallicity. This strong increase of 
aeo at low metallicity is consistent with previous studies, 
in particular those of Local Group dwarf galaxie s which 
model dust emission to co nstrain H2 masses (Lero v et al.J 
I2OIII : IBolatto et aIll201lD . Even though it is difficult to 
parametrize the dependence of aco on metallicity given 
the currently available data the results suggest that CO 
is increasingly difficult to detect at lower metallicities. 
This has direct consequences for the detectability of star- 
forming galaxies at high redshift which presumably have 
on average sub-solar metallicity. 
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